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ABSTRACT 


A  general  principle  for  error  estimation  is  described 
which  can  be  applied  to  different  types  of  partial  differential 
equations.  Particular  attention  is  paid  to  nonlinear  problems. 
With  a  programmed  procedure  based  on  this  estimation  principle, 
error  bounds  are  calculated  for  boundary  value  problems  involving 
the  differential  equation  -Au  +  f(x,y,u)  =  0. 
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1.  The  Problem: 

This  paper  is  concerned  with  a  certain  type  of  differential 
inequalities  which  can  be  used  for  different  purposes.  Here  we  are 
mainly  interested  in  the  application  of  such  inequalities  for  error 
estimations  in  partial  differential  equations. 

For  illustration,  consider  the  boundary  value  problem 


-Au  +  f(x1,x2,u)  =  0 


,  1  2s 

u  =  s(  X  ,  X  ) 


on  G, 

on  r, 


(1.1) 


1 


f 


I 


1  2 

where  G  is  a  bounded  open  domain  of  the  (x  ,x  ) -plane  and  T 
denotes  its  boundary.  Under  appropriate  assumptions  on  f  and  the 


functions 

u ( x1 ,  x2) ,  v( 

x\x2),  it  can  be 

shown 

that 

the  inequalities 

-Au  + 

f(x1,x2,u)  < 

-Av  +  f  (x^,x2,  v) 

on 

G» 

uU^X2)  ^ 

v(x  ,x  ) 

on 

r, 

(1.2) 

imply 

,  1  2s  y 

u(x  ,x  )  < 

v(x  ,x  ) 

on 

G  = 

g ur  . 

We  are  concerned  with  such  implications  for  a  more  general 
problem.  We  consider  a  problem  with  n  unknown  functions  u^(x\ . . . ,xm) 
(i  =  1,2, . .  .,n)  of  ra  independent  variables.  These  functions  shall  be 
defined  and  continuous  on  the  closure  G  of  a  bounded  open  domain  G 
of  the  m-dimensional  Euclidean  space  with  boundary  T.  With  the 


! 


t 

i. 

f 


I 
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notations  x  =  ( x\  . . .  ,xm)  =  ( x*)  ,  u  =  (u\ . . .  ,un)  =  (uk)  =  u(x) 

the  problem  shall  consist  of  p  equations 

F^[u](x)  =  r^(x)  on  B.  ( j  =  1,2, . . .  ,p)  .  (1.3) 

J 

Each  of  these  equations  is  a  differential  equation  of 

at  most  second  order,  given  on  some  subset  B^.  c:  G  .  Besides 

1c 

continuity,  the  functions  u  (x)  shall  have  appropriate  differenti¬ 
ability  properties  depending  on  the  special  problem.  Let  R  denote 
the  linear  set  of  all  u(x)  which  have  these  properties. 

Then,  under  what  conditions  is  the  following  implication  true: 

FJ,[u](x)  ^Fj[v](x)  on  B. 

imply  uk(x)  ^  vk(x)  on  G 

for  u,  v  e  R  . 

Example  1:  The  problem  (1.1)  can  be  written  in  the  form  (1.3) 

using: 

F^[u](x)  =  -Au  +  f(x1,  x2,  u) ,  r^(x)  *  0,  B^  =  G  j 
F2[u](x)  =u(x)  ,  r2(x)  =  s(x1,  x1) ,  B2  =  r  * 

m  =  2,  n  =  l,  p  =  2.  As  R,  one  may  choose  the  set  of  all  functions 
u(x)  which  are  continuous  on  G  and  have  continuous  first  and  second 
derivatives  on  G  . 

Often,  it  is  of  advantage  to  describe  the  problem  in  a  shorter 
abstract  form.  Let  Mu  denote  the  vector 


( j  _  2, .  • .  ,p) 

(k  -  1,2, ... , n) 


(1.4) 


***%  i+imMmm'******** ’*.*&»*+  mm 


r  1 r*,«>-  *  •  **  ntt“+- rtrfcs* tjmw*  ty»»- w*)ffg  tsxxmmu  ft  .iw*e\m.**t.+  mww  -■»*«* ' >  1 


3 

M  u  =  (F1[u],...,Fp[u])  = (F^ [u]) 

which  is  an  element  of  the  set  S  of  vectors  U  =  (U^(x) ,  . . .  ,Up(x) ) 

with  components  U^x)  defined  on  B.  . 

J 

Then,  with  r.=  (rJ),  the  problem  (1.3)  is 
M  u  =  r  $ 

and  the  implication  (1.4)  can  be  written  as 

•  Mu  <  Mv  implies  u  £  v,  (1.5) 

where  inequalities  between  vectors  of  functions  are  defined  in  a 
natural  way,  namely,  as  holding  component  and  pointwise. 

A  problem  (1.3)  satisfying  (1.4)  is  often  called  of  monotonic 
type  [3]  ,  [7],  or  inverse -monotonic  [13].  An  operator  M  with  property 
(1.5)  is  also  said  to  be  inverse -monotonic  (more  precisely:  inverse-isotonic). 
The  reason  for  this  notation  is  that  (1.5)  is  equivalent  to  the  following 
statement.  The  inverse  operator  M  ^  exists  and  is  monotonic  (isotonic): 

U  <  V  implies  M_1U  <  M_1V, 

if  U,V  are  in  the  range  of  M.  Therefore,  an  equation  Mu  =  r  with 
inverse -monotonic  operator  M  has  at  most  one  solution. 

We  will  use  the  following  notations:  Derivatives  with  respect  to  a 
variable  x^  are  denoted  by  subscripts,  for  example,  =  d^u^/d^dx^  . 

If  n=l,  we  will  write  u=u^.  In  this  case,  a 
differential  operator  F^  of  the  type  we  consider  is 
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p^[u]  =F^[u](x)  =F^(x,u,ui,uki)  (i,k,A  =1,2, .  ..,m)  .  (1.6) 

This  equation  is  to  be  understood  in  an  obvious  way*  F*^[u]  may 
depend  on  all  first  and  second  derivatives.  In  (1.6) ,  the  values 
u,ui,ukX  also  k®  considered  as  independent  variables  of  a 
function  F1^ .  For  simplicity,  it  will  always  be  assumed  that  this 
function  F^  is  defined  for  xc  Bj  and  all  values  of  the  other 
variables.  Similar  assumptions  shall  be  made  for  the  other  examples 
of  this  paper. 

2.  Application  to  Error  Estimation 

Suppose  the  given  problem  Mu  =  r  has  a  solution  u**R  and  the 
implication  (1.5)  is  true  for  u,  vcR.  Then,  if  v,  w*R  satisfy 

Mw  ^  r  ^  Mv, 

one  can  conclude  that 

w  ^  u*  <  v. 

Given  an  approximate  solution  cp  with  defect  d[cp]  =  -Mcp  +  r, 
one  can  try  to  get  such  elements  v,w  in  the  form: 
v  =  cp-3z,  w  =  cp+0z,  where  z«R  is  suitably  chosen  and  p  denotes 
a  real  number. 

For  such  v,w,  the  statement  above  is  equivalent  to: 

M(cp-Bz)  -Mcp£  d [cp]  <  M(cp+  gz)  -Mcp  (2.1) 

*  • 

implies 

-0z  <u*  -<p  £  0z. 


ks&rp$sto*  r*wm<*mm**vi**vQ<  ****»«**, 


•»w«ww scmm 
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More  explicitly,  the  inequalities 

Fj[cp-3z]  -FJ’[cp]  <  -Fj[cp] +rJ  £  Fj[cp  +  3z]  -FJfcp]  on  3.  ( j=l,2, . . .  ,p)  (2.2) 

J 

imply  that 

|u*^  -  cp^|  <  3z^  on  G  (k=l,2, . . .  ,m)  (2.3) 

Example  It  For  the  example  of  Section  1,  this  implication  takes 
the  following  form: 

3Az  +  f(x,cp -gz)  -f(x,cp)  ^  Acp-f(x,cp)  <  -3Az  +  f(x,cp+ 3z)  -f(x,cp)  on  G, 

1  o  i  o  (2.4) 

| cp  -  s(x  ,x  )  I  <  3z(x  ,x  )  on  T 

imply 

| ( x1 , x^)  - cp( x"1 , x^)  |  <  3z(x1,x<‘)  on  G  (2.5) 

Error  estimations  of  the  form  (2.2),  (2.3)  are  what  we  are 
interested  in  here.  It  depends  on  the  type  of  the  problem  how  cp  and 
the  error  bounds  can  be  calculated.  Usually,  the  most  difficult  part  is 
to  calculate  an  approximate  solution  cp  with  small  defect.  We  realize 
that  in  practical  applications,  it  might  not  be  the  solution  u*  of  the 
differential  equations  which  is  of  main  interest,  but  other  things,  as 
derivatives,  certain  linear  functionals,  etc.  In  many  cases,  however, 
the  solution  is  of  interest*  and  in  even  other  cases,  bounds  for  the 
solution  may  help  to  get  other  information  which  one  wants  to  have. 


3.  Several  Proofs 


By  different  methods,  several  types  of  problems  have  been  proved 
to  be  inverse -mono tonic.  We  sketch  some  of  the  typical  proofs  without 
stating  all  necessary  assumptions  concerning  differentiability,  etc. 

The  proofs  can  be  applied  to  more  general  problems.  Some  of  the  results 
which  we  will  mention  do  not  exactly  have  the  form  (1.4)  but  involve 
the  <  sign,  also. 

Elliptic  boundary  value  problems:  We  will  prove  (1.2)  under 
different  assumptions. 

1)  The  impllcati on  (1.2)  is_  true  if  on  G: 

f(x,u)  <  f ( x, v)  for  u<v.  (3.1) 

Proof:  Suppose  w  =  v-u  has  a  negative  minimum  at  x=£eG. 

Then  Aw  >0,  and  consequently,  -Au  +  f(x,u)  >  -Av+f(x,v)  at 
x  =  C  in  contradiction  to  the  assumption  in  (1.2). 


In  the  same  way,  one  can  prove  that 
-Au  +  f(x,u)  < -A  v  +  f(x,v)  on  g| 


imply  u  £  v  on 


u  <£  v  on  T I 

if  on  G:  f(x,u)  £  f(x,v)  for  u  £  v. 

Using  the  strong  raximum  principle  one  can  show  that 
2)  The  implication  (1.2)  i_s  true  if  on  G: 


G 


> 
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Proof  [3]  i  [7]:  Let  w  and  K  be  as  before,  and  let  K<-G 
be  any  neighborhood  of  K  such  that  w(x)  ^  0  in  K.  Then,  w 
accepts  its  minimum  at  an  inner  point  of  K  and  satisfies  Aw  >  0 
on  K.  Therefore,  according  to  the  strong  maximum  principle  for 
elliptic  equations,  w(x)  is  constant.  This  is  not  possible 
because  w(5)  <  0,  and  w  >  0  on  T. 

In  another  proof  [2],  a  contradiction  is  derived  by  multiplying 
the  first  inequality  in  (1.2)  by  v,  then  integrating  over  a  small 
enough  sphere  with  center  K  and  applying  Green* s  formula. 

The  statement  (1.2)  can  be  proved  under  even  weaker  assumptions. 

For  example,  in  case  of  a  linear  operator  with  f(x,u)  =q(x)u-g(x) 
variational  methods  have  been  used  [l],  [6].  The  variational  problem 
which  has  -Aw  +  q(x)w=  r(x)  as  its  Euler  equation  is  shown  to  possess 
a  unique  solution  if  on  G 

q(x)  >  -X1  +  6, 

where  6-const  >  0  and  denotes  the  smallest  eigenvalue  of  the 

Dirichlet  eigenvalue  problem  -Au  =  \u  .  Then,  it  is  proved  that  in 
case  r  >  0  a  function  w  which  has  negative  values  does  not  solve 
the  variational  problem. 

The  last  result  can  be  applied  to  the  nonlinear  case  by  using  the 
mean  value  theorem  of  differential  calculus,  assuming  that  df/^u  exists. 

3)  The  implication  (1.2)  is  true  if  on  G: 

|^(x,u)  >  -\1  +  6  . 


(3.2) 
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Another  way  of  extending  the  simple  result  1  is  by  applying 
the  idea  of  the  first  proof  to  w=(v-u)/z  instead  of  w  =  v-u, 
where  z(x)  >0  on  G  .  In  case  of  a  differentiable  function  f, 
one  then  comes  out  with  the  condition 

-Az  4  f^(x,u)z  >0  on  G 

instead  of  (3.1).  Assuming  appropriate  knowledge  of  the  eigenfunction 
corresponding  to  and  choosing  z  close  to  such  an  eigenfunction, 

one  can  get  the  condition  (3.2)  in  this  way,  also. 

Parabolic  problems >  In  order  to  explain  a  typical  proof  for  initial 
value  problems,  we  consider  a  problem  for  one  dependent  variable  u 
( [3]  [4]  [8]  [9]  [10]  [12]  [16])  i 

4)  The  inequalities 

v2 -f(x1,x1,v,v1,vn)  (3.3) 

0  <  x2  £  X; 

0<x2^X,  (3.4) 


0  <  X2  <;  x, 

if  the  function  f  is_  isotonic  (increasing)  with  respect  to  its  last 
variable. 

Proof.  Suppose  that  w  =  v-u  >  0  for  x2  <  C2,  but  w(C1,^2)  =0 


1  2 

u2-f(x  , x  ,u,u1,un)  < 

for  0  <  x^  <  1, 
u  <  v  for  i  x1  =  0,1 ; 

|o  1  x1  <  lj 

together  imply 

u  <  v  for  0  <  x^  <  1, 
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at  a  fixei  point  C  =  (^,£2).  Then,  w2  <  wl  =  WH  -  ® 

at  x=S.  These  inequalities  contradict  (3.3)  because  f  is 
isotonic  in  u^  . 

Hyperbolic  Equati oust  In  hyperbolic  equations,  comparably  few 
results  are  known.  In  some  papers  [4],  [5],  [ll],  [15],  identities 
are  used  which  can  be  derived  by  Green's  formula,  such  as  the  very 
simple  one  (D'Alembert's  formula): 


ff  u12dy1dy2  =  -g[u]  +  u(x)  (3,5) 

B 

x 

where 

y 

g[u]  =|[u(x1,l  -x1)  +  u(l  -x2,x2)]  +  +  U2^y1^1  -y1)^ 

1  -X 

12  2  11 
and  B  consists  of  all  y  =  (y  ,y  )  such  that  1  -y  <  y  £  x  , 

,  1/2.2 
1  -  x  <  y  £  x  . 

5)  The  inequalities 

u]2^  f(x1,x2,u)  <  vl2  4  f(x]‘,x2,v)  (3.6) 

for  l-x2<x1<l,  0  <  x2  <  1, 

for  x2  =  1  -  x^ ,  0  <  x^  <  1  (3.7) 

t  igether  imply 

u  <  V  f_£  1  -  x2  <  x1  <  1,  0  <  X2  <  1  , 

if  the  function  f  is  anti  tonic  with  respect  to  its  last  variable. 


U  S  V 


U1  4  u2  ^  V1  4  v2 
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Proof:  The  inequalities  (3.7)  imply  that  g[u]  <g[v].  Therefore, 
by  integration  of  (3.6)  and  ’ising  the  identity  (3.5)  ,  it  follows  that 


u  + 


// 

B 


f(y1»y2>u)dy1dy2  ^  V  + 


f(y1>y2>v)dy1dy2 


(3.8) 


Suppose  now,  that  w  =  v-  u>0  in  B^,  except  x=  £  where  w(0  =  0. 
Then,  w^  <£  0  and  w^  <[  0  at  x=(.  These  inequalities  contradict  (3.8) 

because  f  is  antitonic  in  u. 


Discussion:  While  some  proofs  use  other  means,  in  most  of  the 

proofs  given  above,  contradictions  are  derived  by  using  conditions  on 
a  function  w  and  its  derivatives  at  a  certain  point  K.  The 
procedures  differ  for  the  different  types  of  problems.  For  the  boundary 
value  problem  (see  1  and  2)  w(5)  is  minimal,  while  for  the  initial 

value  problems,  the  point  K  — in  a  certain  sense — is  the  "first"  point 
where  w(0  =0  (see  4  and  5).  Because  of  this  relation  w(5)  =0,  for 
the  statement  4  no  restrictions  are  required  concerning  the 
dependence  of  the  operator  on  u,  such  as,  for  example,  (3.1).  (Of  course, 
this  pro  of -technical  difference  has  some  deeper  reason,  which  becomes 
apparent  by  the  formula  (3.2)  involving  an  eigenvalue.)  The 
proofs  of  4  and  5  require  some  strong  inequalities  as  (3.4)  ,  (3»7).  Such 
implications,  involving  the  <f  -sign,  usually  can  be  proved  easier  and  under 
weaker  assumptions!  and,  in  fact,  they  do  not  have  the  same  consequences. 

For  example,  the  existence  of  M  \  i.e.  the  uniqueness  of  a  solution  does 


not  follow. 
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The  following  section  contains  a  unifying  approach.  For  all 
types  of  problems  considered  here,  a  weaker  type  of  implication 
involving  the  <-sign  can  be  proved  in  thf  same  way  by  using  conditions 
on  w  at  a  point  K  where  w  is  minimal  and  w=0  (Assumption  I';. 

The  general  type  of  thf  problem,  or  the  properties  of  a  special  problem 
are  then  taken  into  account  by  the  construction  of  a  certain  element  z 
(Assumption  II)  .  Roughly  spoken,  in  this  way  the  assumptions  are  split 
into  "local"  and  "global"  conditions. 


4.  A  Unifying  Approach 


Many  of  the  known  results  about  different  types  of  problems,  and 
new  results  also,  can  be  gained  in  a  unique  way  by  applying  an  abstract 
theorem  on  inverse -monotonic  operators: 


Theorem  [13]  •  Suppose  that  the  following  assumptions  are  satisfied 


I.  For  abritrary  u,  u  e  R, 

u  <  u 
Mu  ^  Mu 

II .  For  a  given  vc  R  there  exists  z  e  R  such  that 


imply  u  ^  u 


(4.1) 


z  >  0,  Mv  ^  M(v  +  \z)  for  X  >  0. 


(4.2) 


Then,  for  all  ueR, 

Mu  <  Mv  implies  u  <  v.  (4.3) 

Remark:  The  theorem  remains  true  if  all  inequalities,  except 


X  >  0,  are  reversed. 
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The  notations  which  occur  in  the  theorem  can  be  defined  in  an 
abstract  way.  For  our  problem,  we  define  an  inequality  u^  v 
between  vectors  of  functions  as  component-  and  point-wise  strict 
inequalities: 

u  <4  v  iff  uk(x)  <  vk(x)  for  all  k  and  x, 


for  which  these  inequalities  have  a  meaning.  Similarly,  Mu^  Mv 
is  defined.  Then  the  assumptions  take  the  following  form: 


Assumption  I:  For  u,  ueR, 
uk(x)  £  uk(x) 

FJ[u](x)  <  Fj[u](x) 


the  inequalities 

(xc  G  j  k“lj  2,  •  •  •  j  n)  y 

(x«B,  j  j  =l,2,...,p) 


(4.4) 

(4.5) 


together  imply 

uk(x)  <  uk(x)  (xcG;  k  =  1,2, . .  .,n)  .  (4*6) 

Assumption  II:  There  exists  z c  R  such  that 

z(x)  >  0  on  G,  FJ'[v](x)  <  FJ'  [  v  +  X  z]  ( x)  (xcB  j=l,  2,  ...,p ;  X  >  0)  (4.7) 

J 

The  element  z  occurring  in  the  error  estimation  (Section  2)  can 
often  be  used  in  Assumption  II,  and  vice  versa. 


We  will  give  an  idea  of  the  proof  for  the  simple  special  case  where 

(1.3)  consists  of  p  =  2  ordinary  equations  for  n  =  2  unknown  numbers 
1  2 

u  ,u  .  (One  may  consider  these  unknowns  as  constant  functions  defined 
on  some  domain  G  )  .  Then,  if  a  vector  v  is  fixed,  the  vectors 
u  ^  v  (u  ^  v)  constitute  a  closed  (respectively  open)  quadrant. 
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First,  notice  that  z  ^  0  because  of  (4.2),  (4.1).  Let  now 
Mu  <  Mv  ,  but  not  u  <  v.  Then,  the  position  jf  the  "points" 
u,  Mu,  V4  Xz  ,  nay  be  as  in  Figure  1.  The  existence  of  a  point 
u  =  v  4- >.  z  such  as  in  this  figure  contradicts  Assumption  I. 


This  proof  shows  that  the  theorem  can  be  generalized  in  different 
ways.  For  example,  (4.1)  need  not  be  required  for  all  u,  u  c  R. 
Moreover,  one  may  consider  a  more  general  curve  v(\)  instead  of 
v  4  X  z  ( see  Figure  l)  . 

.  Local  Assumptions 

The  general  idea  of  proving  Assumption  I  is  as  follows.  Suppose 
(4.6)  is  not  true.  Then,  for  some  k  =  k^,  w^=u^-u^  has  a  minimum 
w  ( K)  =0  at  some  KtS.  This  fact  yields  some  conditions  on  certain 

kQ 

first  and  second  order  derivatives  of  w  at  x-K  .  The  set  of  these 


r 
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conditions,  together  with  the  inequalities  w  (x)  >  0  (4.4),  are 

*i  use  i  to  derive  a  contradiction  to  one  of  the  inequalities  in 


J 


Or 


(4.5).  For  this,  there  must  exist  some  j  .  such  that  F  [u] 

k  U 
K0 

contains  derivatives  of  the  component  u  ,  only.  Moreover,  the 

jQ 

function  F  must  satisfy  some  "monotonicity-condition"  (M-condition)  . 


We  give  some  examples  where  this  proof  can  be  carried  through  if 
the  given  M-conditions  are  satisfied.  It  is  assumed  that  the  occurring 
functions  possess  appropriate  differentiability  properties.  We  will 
not  specify  the  set  R  of  those  functions  in  each  case.  The  examples 
can  be  generalized,  in  particular  with  respect  to  the  boundary  conditions. 


Example  2:  Generalization  of  Example  1,  boundary  value  problem  for 
one  unknown  function  u  : 

F1[u](x)  =  r1  ( x)  for  x  e  G, 

F2[u](x)  =  r2(x)  for  x*r 

with 

F1[u](x)  =  F1(x,  u,  u^  ukj£)  (i,k,X=l,2,...,m)  (5.1) 

F2[u]  (x)  =  u.  (5.2) 


M-condition: 


F1(x,u,ui,ukx)  >  F1(x,u,ui,uki  +  wkx) 
for  xeG,  all  and  any  m  x  m-matrix 

("kW  ^  °> 

i.e.  for  any  symmetric  positive  semi -definite  matrix  (wkx)  • 


(5.3) 
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Example  2a:  An  important  special  cas<  jf  (5.1)  is  the  quasi - 

linear  operator 

n 

*u]  -  -  V  akx(x,u,ui)ukj.  4  e(x,u,u.) 
k,*=l 

with  a  coefficient  matri*'  satisfying  the  M-c  or.ditiun: 

(akx(x,i;,ui))  >  0. 

In  xe  G,  and  all  u,u^  • 

For  example,  this  c  ndition  is  satisfied  for  the  equation 
of  minimal  surfaces: 

F1[u]  =  -[1  -f  (ux)  2]un  i  2u]_u2ui2  -  rl  -»  (n2)  2]u22  =  0, 

2  9  2 

and,  if  (u^j  +  (u2)  <  c  >  ^or  equation 

F1[u]  =  -[c2  -  (u1)  2]un  +  2j1u2u12  ~  '  "  (i:2)  21u22  =  °> 

describing  the  two-dimensional  steady,  irrotational  fl  iw  of  a  compressible 
fluid. 

Example  3  ’•  Initial-boundary  value  problem  for  one  unknown  function: 


R 

* — i 

n 

u 

'R 

r~k 

for  x  e  1  U  Tp  , 

(5.4) 

F2[u](x)  =r2(x) 

for  xeE, 

with  F2  as  in  (  5.  2)  and  F^  a  special  case  of  (5.1): 

F1  [u]  =uk -f(x,u,ui,ukx)  (i,k,X  =  l,2,...,m-l)  .  (5.5) 

G  shall  be  a  cylinder  of  points  x,  such  that  ( x^ , . . . ,xra~^)  is 
element  of  an  open  bounded  domain  G,  and  0  <  xm  <  X;  is  the  set 

of  boundary  points  with  (x\...,xm  e  G,  xm  =  X,  and  =  f  -  r2  • 
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M-Condltlon:  Condition  (5.3)  applied  to  ( 5.5)  . 

For  example,  the  equation  of  heat  conduction 

aU4  ~  "  ^2^  2  ‘  (^3)3  +  f  =  0 


with  positive  thermal  properties  a(x,u),  k(x,u),  and  heat  genera ' 1  on 

f(x,u)  can  be  written  as  (5.4)  such  that  the  M-condition  is 
satisfied. 

Example  4«  Semilinear  hyperbolic  systems  of  first  order  in  normal 


form: 


n 

Fj[u](x)  =  -V  orJi 
k=l 


(x)u^  +  fJ(x,  u) 


=  0  for 


x  c  g  u  r . , 

J 


(5.6) 


Fn+^[u](x)  =  uJ'(x)  =sJ'(x)  for  xtrn+j  =  T  -  T 


where  T-  denotes  some  part  of  I\ 
J 

M-Condition: 


i)  f^XjU1,.  .,,un)  >  f^(x,v^,...,vn) 


uK  <  vk  with  k  ^  j 


—  i  i 

1  uJ  =  vJ 


(xe  G  UT.j  j  -1,2,  ...,n) 

J 


ii) 


For  each 


x,rr 


the  ” directional  derivative” 


> 


is  a  derivative  into  the  closed  domain  G. 


Example  4a  ( Wave  Operator)  1  The  differential  equation 
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V12  -1  f(x1,X2,V,T1,V2)  =  0 

can  be  transformed  into  a  system  (5.6) 

Jr  1  1,1  2  3 

r  [u]  =  +  u2  -  u  -  u  =0, 

»-,2r  i  2  ,  ,,/  1  2  1  2  3\  ^ 

F  [u]  =  U2  4  f(x  ,  X  ,u  ,u  ,u  )  =  0, 

-3r  1  3,-/1  21  2  3\ 

F  [u]  =  u ^  +  f(x  ,x  ,u  ,u  ,u  )  =  0. 

The  first  part  of  the  M-condition  is  satisfied  if  f  is  antitonic 

(decreasing)  with  respect  to  the  last  three  variables.  For  the  usual 

boundary -initial  value  problems,  the  second  part  is  satisfied  also.  A 

simple  example  is  the  characteristic  initial  value  uroblem  for  the  domain 

0  <  x\x2  <  1  with  u1  given  for  x^  =  0  and  for  x2  =  0( Tj)  >  u2 
13  2 

given  for  x  =  O(T^)  ,  u  given  for  x  =  O(r^)  .  (The  result  5  in 
Section  3  can  be  obtained  by  applying  the  abstract  theorem  to  the  integral 
operator  occurring  in  (3.8)). 

Other  Examples:  The  examples  given  above  can  be  generalized.  Moreover, 
it  is  often  possible  to  transform  a  given  problem  such  that  Assumption  I 
becomes  satisfied.  We  describe  some  of  the  involved  ideas  using  examples. 


1) 


Consider  a  problem  (1.3)  where,  for  x  «  G,  the  equations 


F1[u]  =  (u1  +  u2)u^  +  u^  -  u2  =  r*(x) 

_,2r  I  2  1  2/  \ 

F  [u]  =  u2  -  u  =  r  (x) , 


(5.7) 


are  given.  These  equations  do  not  have  the  form  necessary  to  derive  a 

contradiction  to  w^(0  =0  at  a  point  C  e  G,  because  F^  is  not 

.  .  .  2  12 
antitonic  with  respect  to  u  .  Therefore,  introduce  new  variables  U  ,U 


and  replace  (5.7)  by  the  system 


i 
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J.lr  h  11  ,,2/  1N+  2/  lx-  1  ,  2  _  1/  x 

r  [u,l]  "  u  u.^  -  U  (u^)  -  u  (u^)  +  u2  -  u  -  r  (x), 

F12[u,U]  =  -Fnf-l,-u]  =  -r1  ( x) , 

F21[u,U]  =  F2[u]  =  r2(x),  F22[u,U]  =  -F2[-U]  =  -r2(x), 

where  2f+  =  |f|  +  f,  2f~  =  |f|  -  f.  For  x  c  G,  this  system  has 
all  necessary  properties,  and,  for  =  -u^,  it  is  equivalent  to 

(5.7).  As  approximations,  respectively  bounds,  for  the  new  variables 

K  k  k  k 

U  ,  one  can  use  -cp  ,  respectively  -<p  ±  8z  .  Then,  the  correspond¬ 
ing  inequalities  (2.2)  consist  of  equivalent  pairs,  which  do  not 
explicitly  involve  the  positive  part  f  ,  or  negative  part  f  of  any 
function,  but  |f|  instead, 

2)  If  a  system  (1.3)  does  not  have  "normal  form",  i.e.  if  not  all 
equations  contain  derivatives  of  one  variable  only,  one  may  get  this 
normal  form  by  transformation  of  the  dependent  variables. 


Consider  a  quasi -linear  system  in  matrix  notation: 

A(x,u)u^  +  u2  =  0  on  G,  together  with  appropriate  boundary  conditions. 

Suppose  there  exists  a  nonsingular  matrix  ®(x,u)  and  a  nonsingular 

diagonal  matrix  D(x,u)  such  that  C  ^A®  =  D.  Then,  the  given  system  is 

equivalent  to  D*>  +  '■*>  \i2  =  0.  One  can  introduce  new  dependent  variables 

k  -1  -1 

v  such  that  <*>  u^  =  |v^  +  c,  $  u2  =  fv2  +  d,  if  the  function  i|i(x,u), 
and  the  vectors  c(x,u) ,  d(x,u)  can  be  determined  such  that  (v^)2=  (v2)  . 
Then,  the  transformed  system  D(  fv^  +  c)  +  ( fv2  +  d)  =  0  has  normal  form. 

For  example,  use  v  =  <®~^u  if  G>  does  not  depend  on  u. 
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In  the  case  of  the  equations 


2  1  12,  1  _  n 

u  +  u  =  0, 


2  1,  122,  1  2  .  A  ...  2_  r  /  1\  -i  2 

c  u^  +  u  u  u^  +  u  u^  =  0  with  c  =  [c(u  )j 

which  describe  the  one -dimensional  (i  sen  tropic)  flow  of  a  compressible 
fluid,  a  system  in  normal  form  can  be  gained  in  the  indicated  way  for 

the  variables  v^  =F(u^)  +  u^,  v^=-F(ux)  +  u^,  where  F(p)  =  do. 

* 


6.  Constructing  z 

Assumption  II  is  not  very  restrictive,  in  the  sense  that  something 
"not  very  much  weaker"  must  be  required.  In  fact,  for  linear  problems, 
Assumption  II  is  necessary,  in  general.  For,  if  a  linear  operator  M 
is  inverse -monotonic,  and  if  for  some  r  y  0  the  equation  Mu  =  r  has 
a  solution,  then  this  solution  u  =  z  satisfies  ^4.2) 


The  element  z  can  be  constructed  for  certain  large  classes  of 
problems. 

Example  2:  If  the  function  F^  in  (5.1)  is  strictly  isotonic  with 

respect  to  u,  then  z  =  1  is  appropriate.  More  complicated  functions 
z  yield  weaker  restrictions  on  F^.  For  example,  the  function 


/ 


-P(r0-t) 

z=  fte  dt,  has  been  used  [13],  where  r^  is  the  radius  of 


an  open  sphere  K  3  G  and  r  denotes  the  Euclidean  distance  of  its 


center  from  x. 
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Example s  3  and  4a?  For  initial  value  problems  often  very 

simple  functions  z  are  satisfactory.  In  Example  3>  Assumption  II 

m 

is  satisfied  with  z  =  e“X  ,  I!  large  enough,  if  f  obeys  a  certain 
one-sided  Lipschitz-condition.  In  Example  4a,  a  similar  condition 

is  sufficient  if  zk  =  e^ X  +X^  is  used  (k  =  l,2,3). 

In  properly  exploiting  the  theorem  one  can  also  gain  conditions 
of  more  theoretical  nature,  like  (3.2).  Consider,  for  example,  the 
problem  (1.1)  withf  =  q(x)u  -  g(x),  and  assume  that  T  and  q(x)  are 
sufficiently  smooth.  Then,  for  any  such  f  satisfying  (3.^)»  and 

0(x)  el,  s(x)  el,  the  problem  (1.1)  has  a  solution  z.  For  example  z 
for  q  e  1.  Thus,  for  q  =  1,  the  Assumption  II  is  satisfied.  But  it 
is  also  satisfied  for  any  other  f=qu  satisfying  (3.2),  because 
q(x)  can  be  connected  with  q  si  by  a  curve 

q(x,t)  =  tq(x)  +  (1-t)  (0  t  ^  l) .  The  corresponding  z(x,t)  depend 

continuously  on  t,  and  z(x,0)  >  0  on  G.  If  z(x,l)  >  0  on  G  would 
not  be  satisfied,  then  for  some  te(0,l)]i  z(x,t)  ]>  0  on  G,  but  not 
z(x,t)  >  0  on  G.  This  contradicts  Assumption  I.  Obviously,  this  is 
a  very  special  case  of  far  more  general  results. 

7.  A  Program 

The  preceding  sections  have  shown  that,  in  principle,  the  method  of 
error  estimation  of  Section  2  can  be  applied  to  many  types  of  problems. 

Of  course,  for  a  concrete  problem,  usually  a  lot  of  additional  considera¬ 
tions  are  necessary  to  make  the  method  work.  To  investigate  the  practical 
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application,  a  program  has  been  written  for  the  problem  (1.1).  The 
general  idea  is  as  fol  ows.  (For  convenience,  we  use  partly  different 
notations  than  before.  For  example,  we  write  (x,y)  instead  of 
(x  ,x  );  and  subscripts  do  no  longer  denote  derivatives.) 


Step  D i  Calculating  Approximations  by  the  Difference  Method. 

(This  step  is  not  needed  for  linear  problems.)  The  ordinary  difference 
method  is  applied  to  calculate  approximate  values  u^  at  the  net-points 
(x^,y^)eG  of  a  rectangular  net.  If,  for  example,  the  condition  (3.2) 
is  satisfied  and  the  mesh  width  is  not  too  large,  the  nonlinear  difference 
equations  can  be  solved  by  the  iterative  procedure  of  Picard  or  I.ewton. 

We  have  restricted  ourselves  to  cases  where  the  Picard  procedure  converges 
(|fj  <^-6)  ,  arid  we  have  solved  the  nonlinear  systems  by  a  combination  of 
Picard's  procedure  and  the  point-overrelaxation  method. 


Step  A:  Calculating  an  Approxlmati on  cp(x,y)  .  A  development 
cp  =  cp^  +  +  •••  +  ayp^  with  properly  chosen  functions  cp^(x,y)  is 

set  up.  For  calculating  the  constants  ,  the  defect 
d^[cp]  =  Acp  -  f(x,y,cp)  is  replaced  in  the  net-points  (x.  ,y.)  by  the 
linear  approximation 


d1[cp]  =  i>cp  -  f(x,y,u)  -  fu(x,y,u)  (cp-u)  with  u^,^.)  =  u^  . 

In  order  to  get  "small"  defects  d^[<r]  on  G,  and  d^[ <r]  =  -cn  +  s  on  r,  an 
orthogonalization  method  is  applied} 


f  ^  2  ~  0  (k  =  l,2,...,m), 


(7.1) 
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where  ^  and  0^  are  properly  chosen  functions,  the  constants 

and  6^  are  only  introduced  for  computational  convenience,  and 

(  |  )  ^ ,  (  *  )2  deno^e  discrete  inner  products.  For  example, 

(u,v),  =  l  w.  .u(x.,y .)  v(x.,y.)  with  given  weights,  involving  all  net- 
^  J  ^  J  *  J 

points  in  G.  The  linear  system  (1.7)  for  the  ar^  is  solved  by  an 
elimination  method. 

Step  E t  Calculating  an  Error  Bound.  For  a  chosen  function  z,  one 
determines  a  constant  3q  such  that 

I  <  P0(-Az  +  fu(x,  y,  <p)z}  on  G,  ld2[<P]  I  <£  PQz  on  T. 

Practically,  this  is  done  for  points  (x,y)  in  a  finer  net.  Then,  if 
3q  is  small  enough,  the  desired  inequalities  (2.4)  hold  for  a 
number  3,  somewhat  larger  than  3Q,  say  3  =  1.013 }.  This  has  to  be 
checked.  If  this  is  so,  then  the  error  estimation  (2.5)  holds. 
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8.  numerical  Examples 


The  pz’ogram  of  Section  7  has  been  used  to  calculate  approximations 
and  corresponding  error  bounds  for  the  following  problems: 


1. 


2. 


3. 


4. 


-Au  =  1  for  |x|,  |y  |  <  1, 
U  =  — ~r(  COS  TTX  +  COS  TTy) 

2ri 


on  the  boundary. 


-Au  =  1  for  |x| , |y ]  <  1, 
u  =  0  on  the  boundary. 

-Au  =  eU  for  |x| , |y |  <  1, 
u  -  0  on  the  boundary. 

-Au  =  -eU  for  |x|, |y|  <  1, 
u  =  0  on  the  boundary. 


The  problems  1  and  2  were  mainly  calculrtea  to  check  how  the  method 
works,  before  starting  the  nonlinear  problems  3  and  4.  The  general 
program  is  constructed  to  handle  more  complicated  problems.  However,  we 
did  not  accomplish  to  compute  more  examples  during  the  available  time. 

In  all  of  these  problems  we  chose  a  square  net  with  mesh  width  h  =  0.04. 

In  Problem  1,  we  used  the  development 

cp  =  <p0  +  (1  -  x2)(l  -  y2)^.^  +•**+  amum)  (8.1) 


with  polynomials  having  appropriate  symmetry  properties: 
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,  2  2  2  2  (  q  n\ 

1,  x  +  y  ,  xy  (0.2) 

and 

(p  =  — “t( cos  ttx  4  cos  ny)  . 

J  2n2 

As  functions  there  were  used  orthogonal  polynomials 

P0,  P2(x)  +  P2(y) ,  P2(x)P2(y),...  (8.3) 

where  P^(x)  is  proportional  to  the  l  Legendre  polynomial.  Finally 
~  1>  6^-0  (k  =  1,2, . . .  ,m)  ;  w^  ^  —  1  ( i >  j  —  1  >  2,  • .  •  ,ro) ,  (8.4) 

and  in  Step  E 

z  =  2  -  (x2  +  y2)  .  (8.5) 

In  this  way,  we  got  the  following  error  estimations  for  Problem  1. 

For 

m  =  4:  |u*(x,y)  -  <p(x,y)  |  <  0.0093  (2  -  x2  -  y2)  (  |x|,  |y|  <  l) 

where  cp(0,0)  =  0.326  738; 

m  =  6;  !u*( x,y)  -  <p(x,y)  |  <  0.0064  (2  -  x2  -  y2)  (  |x|,  |y|  <£  l) 

where  cp(0,0)  =  0.326  834; 

m=9:  |u*(x,y)  -  cp(x,y)  |  <  0.000270  (2  -  x2  -  y2)  (|x|,|y|^l). 
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In  a  first  run  we  had  tried  with  the  polynomials  in  (8.2)  as 
functions  ^  instead  of  (8.3).  But  then,  the  linear  systems  which 
determine  the  constants  became  ill-conditioned  for  m  =  6  and 

m  =  9»  and  the  corresponding  approximations  cp  were  less  accurate 
than  the  approximation  for  m  =  4. 

For  Problem  2,  we  did  not  apply  the  approximation  and  estimation 
procedure  immediately,  because  Au  is  discontinuous  at  the  corners 
of  the  domain.  The  singularities  were  removed  by  introducing  a  new 
variable 

v  =  u  +  p,  (8.6) 

where  p(x,y)  satisfies 

Ap  =  0  for  all  x,y,  except  at  the  corners, 

Ap  =  1  at  the  corners. 

The  function  p(x,y)  consists  of  four  summands  of  the  type 
-1  2 

tt  Im(z  log  z) ,  each  of  them  belonging  to  one  of  the  corners.  For 
example,  z=l+x+i(l+y)  for  the  corner  (x,y)  =  (-1,-1) .  More 
explicitly, 

P(x,y)  =  ^(Px  +  P2  +  P3  +  P4) 

with 

P2  =  q(x,y),  p2  =  q(-y,x) ,  p^  =  q(-x,-y),  p4  =  q(y,-x) 

and 

q(x,y)  =  (1  +  x)(l  +  y)  log  [(1  +  x)2  +  (1  +  y)  2] 

+  [(1  +  x)2  -  (1  +  y)2]  arctg  . 


26 


The  transformed  Problem  21  then  is 

-Av  =  2  for  |x| , |y |  <  1, 

v  =  p(x,y)  at  the  boundary.  (8.7) 

The  procedure  in  Section  7  was  applied  to  this  problem.  We  again  used 
a  development  (8.1)  with  as  in  (8.2),  but 

cp3  =  £[H(x)  +  H(y)  -  H(l)  ] 

with 

H(x)  =  h(x)  +  h(-x) 

and 

h(x)  =  2(1  +  x)  log  [4  +  (1  +  x)2] 

+  [4  -  (1  +  x)  2]  arctg  £(1  +  x) 

This  function  9^  satisfies  the  boundary  condition  (8.7) . 

With  functions  ^  as  in  (8.3),  the  quantities  in  (8.4),  and  z 
in  (8.5),  we  got  the  following  error  estimation  for  m  =  4s 

|v*(x,y)  -  <p(x,y)  1  =  |u*(x,y)  +  p(x,y)  -  q>(x,y)  |  < 

<  0.000  002  (2  -  x2  -  y2)  (|x|,,y|^l) 

where 

9(0,0)  =  1,177  227  9, 

p(0,0)  =  0.882  542  4, 

9(0,0)  - p( 0, 0)  =  0.294  685  5. 
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Of  course,  the  error  bound  does  not  take  into  account  the  rounding 
errors  which  occurred  during  the  calculation  of  cp(x,y)  according 
to  (8.1). 

For  the  linear  Problems  1  and  2,  we  could  also  have  used  functions 
cp  which  satisfy  the  differential  equation  instead  of  the  boundary 
condition.  (Then,  =  0,  6^  =  1  (k = 1,2, . . . ,m) ) .  3ut,  this  is  not 

possible  for  the  nonlinear  problems  3  and  4>  and,  as  mentioned  above, 
the  linear  problems  mainly  served  to  check  the  procedure. 

The  Problem  3  also  was  transformed  using  (8.6)  before  applying 
the  approximation  and  estimation  procedure.  The  transformed  Problem  3*  is: 

-LV  -  e~p^x,y^eV  for  |x|,|y|  <1, 
v  =  p(x,y)  on  the  boundary. 

In  Step  A  and  Step  E,  we  used  the  same  quantities  as  in  Problem  2’ . 

In  this  way,  we  got  the  following  results  for  m  =  6: 

|v*(x,y)  -cp(x,y)  |  =  |u*(x,y)  tp(x,y)  -cp(x,y)  |  <  0.001  02(2  -  x2  -  y2) 
for  |x|,  |y |  <  1, 

where:  cp(0,0)  =  1.278  072  3, 

9(0,0)  -  p(0,0)  =  0.39  5  529  9, 

and  the  approximations  v(0,0)  for  v*(0,C), 
u(0,0)  for  u*(0,0),  obtained  by  the  difference  method,  are 
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In  Problem 


was  used. 

|v*(x,y) 

where: 


v(  0, 0)  =  1 . 277  974  7 , 

u(0,0)  =  v(0,0)  -  p(0,0)  =  0.395  432  3. 

4  we  proceeded  similarly,  except  that  now  the  transformation 
v  =  u  -  p 

The  results  for  m  =  6  were: 

-  ?U,y)!  =  !u*(x,y)  -  p(x,y)  -  cp(x,y)|  <  0.00038  (2  -  x2  -  y2) 
<p(0,0)  =  -1.127  658  9, 

cp(0,0)  +  p(0,0)  =  -0.245  116  5, 

v(0,0)  =  -1.127  576  2, 

u( 0,0)  =  v(0,0)  +  p(0,0)  =  -0.245  033  8. 
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